Noise-resistant optimal spin squeezing via quantum control 
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Entangled atomic states, such as spin squeezed states, represent a promising resource for a new 
generation of quantum sensors and atomic clocks. We demonstrate that optimal control techniques 
can be used to substantially enhance the degree of spin squeezing in strongly interacting many-body 
systems, even in the presence of noise and imperfections. Specifically, we present a time-optimal 
protocol that yields more than two orders of magnitude improvement with respect to conventional 
adiabatic preparation. Potential experimental implementations are discussed. 

PACS numbers: 



Quantum squeezed states are among the most interest- 
ing examples of entangled states. In quantum metrology 
they allow for measurements with an improved precision, 
ultimately limited only by the Heisenberg limit. Sinee the 
early theoreticalproposals to realize them with non linear 



early tneoreticalproposals to realize them with non linear 
interactions [2|, |3[ , spin squeezed states have been imple- 
mented in several experiments. Specific examples include 



generation of spin squeezed states in cavity QED 

in trapped ions through shared motional modes 

or using non linear interferometry with a Bose-Einstein 

condensate [9|. 

In this Letter we demonstrate that optimal control can 
be effectively employed to produce highly squeezed spin 
states in many-body quantum systems, drastically reduc- 
ing the impact of relaxation and decoherence. We em- 
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ploy the CRAB technique 



11 



to optimally control 



the evolution of a collection of N two-level systems mu- 
tually coupled through a non linear (i.e. quadratic) inter- 
action. We calculate optimized evolutions occurring on 
time scales several orders of magnitude shorter than the 
corresponding adiabatic evolutions, with a speed-up in- 
creasing with the system size. Such a speed-up translates 
directly into an enhanced robustness of the squeezing in 
the presence of noise, as schematically depicted in Fig. [T] 
We illustrate this enhanced robustness by modelling two 
practical experimental implementations of squeezed state 
preparations: cavity QED and trapped ions I6|,l8|. 

One approach to spin squeezed states is based on the 
so called one-axis twisting protocol, consisting in letting 
evolve a collection of two-level systems under the effect 
of a collective non linear interaction |3| , described by a 
Hamiltonian of the form 



HsM = w J^ + xJ^ 



(1) 



for the global spin variable J (see below). Here to is 
precession frequency and x is a strength of nonlinear in- 
teraction. The relative simplicity of the one-axis twist- 
ing scheme has been at the basis of its ubiquitous pres- 
ence in squeezing experiments; however such a scheme 
is known to be non optimal [3|, the spherical nature of 
the angular momentum phase space limiting the maxi- 




FIG. 1: Upper panel: initial state (left) and final highly 
squeezed state (right) for a system of A'^ = 100 spins. Lower 
panel: adiabatic (red) and optimal (black) driving fields x 
generating the maximally squeezed state shown above; the 
effect of the noise (big blue arrow) increases with the total 
evolution time. 



mal squeezing achievable. Such a bound is intrinsic for 
the one-axis twisting protocol with fixed x- S0rensen and 
M0lmer proposed a protocol based on adiabatic evolution 
to steer a system into maximally squeezed states [12J |. 
This procedure has been implemented experimentally in 
low-dimensional systems, see for instance Ref. [13|. Un- 
fortunately the required evolution time, which is propor- 
tional to the inverse square of the minimum spectral gap 



A encountered during the evolution. Tad oc A~^ [14l |. 
scales unfavourably with a system size. This makes adi- 
abatic evolution significantly exposed to external noise: 
typically in many-body systems the gap closes with in- 
creasing system size N, which implies a dramatic in- 
crease of the time required for adiabatic evolutions for 



large N. Previous studies have demonstrated that opti- 
mal control is a powerful tool to drastically reduce the 
time needed to perform a many-body quantum evolu- 
tion [ll|, |l5[ . In particular the Chopped RAndom Basis 
(CRAB) technique offers an efficient way to implement 
optimal control, based on an expansion of the control 
field onto a truncated basis 10|, lll[ . Recently it has been 
shown that optimal control allows for reaching the Quan- 
tum Speed Limit (QSL), the minimal time required by 
physical constraints to perform a given transformation, 
in spin chains [15|, |l6[ , cold atoms in optical lattices [l7| , 
Bose-Einstein condensates in atom chip experiments and 
in crossing of quantum phase transitions [l9|. Indeed, 
CRAB control makes it possible to reduce the time of the 
transformation down to the QSL, which scales as 1/A, 
obtaining a quadratic speedup of the protocol with re- 
spect of the adiabatic one. In this work we show that 
this method is successful also in drastically reducing the 
preparation time for maximally spin squeezed states, as 
illustrated in Fig. [H thereby significantly enhancing the 
process' robustness to reahstic noise sources. 

Model — A collection of N two-level atoms having 
(pseudo)spin Si can be described in terms of the global 
spin variable J = J2i=i^i' '^ith \J\ — N/2 and z- 
component Jz representing the population imbalance be- 
tween the two atomic internal states. In Ramsey spec- 
troscopy experiments, the measured signal M yields the 
mean global angular momentum pointing along the z- 
axis, M = (Jz), while the noise is given by the un- 
certainty in one of the orthogonal components AJi ~ 
\/{Jf) — {Ji}'^, i = x,y. In spin squeezed states, the lat- 
ter is below the standard quantum limit, i.e. AJi < 
\{Jj)\/2 for i ^ j G {x,y, z}. The squeezing parameter S, 
is defined through the signal to noise ratio as 



(2) 



2JAJ^ 



Squeezed states satisfy the condition ^ < 1, which im- 
plies entanglement in the system. The ideal states for 
spectroscopy experiments are those minimizing AJ^ for 
sufficiently large values of the signal, i.e. M ex N. The 
problem of finding the optimal squeezed state can be re- 
cast into the search for the ground state \ipoix-: ^)) of the 
Hamiltonian Eq. [U where lo is constant and negative and 
the non-linear interaction xit) is now taken to be tunable 
in time (26| . (From now on we set h = 1 and time is 
measured in units of l/|a;|.) Adiabatic evolution under 
HsM automatically produces optimal squeezed states, as 
follows: At the time i = one takes x(0) = and the 
system is prepared in its initial ground state \ijJo{0,N)), 
the coherent state \Jz ~ J) with ^ = 1. Then adiabat- 
ically increasing x(Oi the system evolves following the 
instantaneous ground state |V'o(x(*)7-^)) oi Hsm, yield- 
ing exactly the family of states with optimal squeezing 
at a given value of M (see Fig. [IJ . 
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FIG. 2: Scaling with the size of the total evolution time T for 
the adiabatic (lad ~ 7 ■ 10~^, red triangles) and the optimized 
dynamics {lopt = 5 ■ 10^'*, black circles). Numerical fits for 
30 < iV < 150 (dashed lines) result in Tad oc iV^'^^ and 
Tapt oc iV°■^^ Inset: Scaling of ^^ a fit gives ^^ ~ 2.1/N^-^\ 



Optimization in the absence of noise — Wc first inves- 
tigate the properties of the Hamiltonian Hsm in Eq. ([T]) 
to identify target squeezed states that can be reached 
via adiabatic evolution. We calculate the time required 
to perform an adiabatic transformation from the initial 
state into the target and its scaling with the system size 
A^. Subsequently, we apply optimal control to determine 
the dynamics (neglecting for the moment decoherence ef- 
fects) leading to the same target state in a much shorter 
time. Finally, we compare the optimized evolution with 
the adiabatic one. 

As previously mentioned, squeezed states suitable for 
quantum metrology should have sufficiently strong sig- 
nal M . To fulfill this requirement we choose (throughout 
the whole work) M = J/V2 = 0.707 J, i.e. M (x N. 
Then we find the value Xm{^) of the interaction such 
that \ipo{xM,N)) has (Jz) = M for a given N. The 
inset of Fig. ([2]) shows the corresponding value of the 
ground-state squeezing for varying N: a power-law fit 
^2 = A/N^ for 30 < A < 150 yields A = 2.1 ± 0.05 and 
B = 0.94 ± 0.01, compatible with the Heisenberg limit 
^2 ex: N^^. This means that we have identified a class of 
states \tpo{xMT^)) '^ith the desired characteristics. We 
can now take those states as a target for the optimiza- 
tion, to achieve constant intensity of the signal M and 
maximal squeezing ^ for any given system size N. 

As discussed above, the system is initially prepared in 
the coherent state |^o(0,A^)) where all spins are polar- 
ized along the positive ^-direction and ^^ = 1, and we 
aim at reaching the goal state Itpo) = |V'o(Xmj-^)) ^^~ 
ter an evolution time T. The initial and target state for 
the case N = 100 are depicted in Fig. [T] (upper pan- 
els). For the adiabatic case, evolution is computed using 
a linear ramp xit) = Xiyit/T. Comparing the resulting 
final state |'0(T)) with the goal state yields the infidelity 
/ = 1-\{iPo{xm,N)\tJj{T))\'^. Fig.Hshows, as a function 
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FIG. 3: Final squeezing ^^ as a function of the size A^, for 
V = 500 for the adiabatic (red triangles) and optimal (black 
circles) dynamics, subject to random telegraph noise with am- 
plitude Ka ~ Kp — 0.05 (empty symbols) and Kc, — Kp — 
(full symbols). Data have been averaged over 24 instances of 
disorder. 



of the size N, the time Tad needed to reach a given infi- 
dehty value lad via adiabatic evolution (red triangles) . A 
fit T = AN^ for 30 < iV < 150 gives A = 0.31 ±0.01 and 
S = 1.95 ± 0.01, in agreement with the prediction of the 
adiabatic theorem Tad ~ 1/A^ ^ 1/N^. We then apply 
the quantum optimal control CRAB algorithm [10|, lll| to 
find the time Tgpt needed by an optimal transformation to 
reach an infidelity lopt ■ More precisely we write the driv- 
ing field in the form x{t) = Xm[^+K^) ^]Li o.j sin(a;jf)+ 
hj cos{u)jt)\t/T ^ where X{t) ensures constant boundary 
conditions, Uj = 2tt/T{1 + rj), tj is a random number, 
and Hf ~ 0(10), and we look for the optimal correction 
(i.e. the coefficients a, b) such that the infidelity is min- 
imised for a given time (for details on the algorithm and 
of its complexity see [llj, l20[). A typical result is shown 
in Fig. [T] (lower panel), while the scaling of the optimized 
evolution time Topt as a function of the size N is shown 
in Fig. [2] (black circles). A power-law fit Topt = AN^ for 
30 < TV < 150 gives A = 0.06 ±0.01 and B = 0.93 ±0.04, 
consistent with our conjectur about the QSL (see above). 
This shows that optimal squeezing preparation results in 
a quadratic improvement in the scaling of the preparation 
time as a function of the system size, while additionally 
reducing the total evolution time by at least two orders 
of magnitude. 

Our discussion up to this point neglected completely 
the effect of noise, which of course is a major concern 
in a real experiment. Therefore, in order to test the ro- 
bustness of the protocol, we simulate the dynamics of 
the system in the presence of noise. We will consider two 
noise models, as different experimental implementations 
of squeezed spin states are affected by different kinds of 
noise. Wc will show that optimized protocols work also 



in the presence of these types of noise, and that they are 
much more resilient to noise than adiabatic protocols. 

Effect of classical noise — A typical situation in which 
classical fluctuations of an external field occur, reflect- 
ing in random fluctuations of the interaction strength, 
is found in trapped ions, also relevant for metrological 
applications [8|. In trapped ion systems, a global ran- 
dom magnetic noise is expected to be the most relevant 
source of disturbance 18| . We include it in our simula- 
tions by adding random classical telegraph noise to the 
control field. We then study the evolution induced by 
the Hamiltonian 



H = xit)[l + Ko,a{t)]Jl + u[l + K^mUz 



(3) 



where a(t),/3(t) are random functions of the time with 
a flat distribution in [—1, 1], changing random value on 
average with frequency v. The case la = Ip = corre- 
sponds to a noiseless evolution of Eq.([T|). In Fig. [3] we 
compare the effect of the noise on the final squeezing ob- 
tained by varying x(i) either linearly in time (empty red 
triangles) or according to the optimized protocol (empty 
black circles). The squeezing ^^ is plotted as a function 
of the size TV, for i' = 500 and for an intensity of the noise 
Ka = Kp = 0.05. As shown in Fig. [Sj the noise effect 
is stronger for larger system sizes, very quickly destroy- 
ing the squeezing for the slow linear (adiabatic) protocol. 
The reason is simple: as shown in Fig. [2l for large sizes, 
e.g. TV > 100, the adiabatic evolution time is three or- 
ders of magnitude larger than the optimized one. Vice 
versa the fast optimal driving turns out to be robust even 
at large sizes and relatively high intensities of the noise, 
resulting in a final squeezing almost equivalent to that 
obtained via the adiabatic process in the absence of noise 
(full red triangles). 

Effect of quantum noise — Finally we discuss a noise 
model suitable for the description of QED experiments 
[61 , in which the effect of the noise is treated through 
the formalism of the master equation. In cavity QED, 
relaxation of the atomic levels towards the ground state 
and leakage of photons outside the cavity are the most 
relevant source of dissipation [J|. In order to estimate 
the effect of the noise in a realistic system, we derive 
the Hamiltonian of Eq. ([T|) from a microscopical model. 
We consider a collection of TV three level atoms with two 
stable ground states \a) and \b) and an excited state |e), 
in an optical cavity; the ground state energy splitting 
is given by Uab and the relevant cavity mode has a fre- 
quency LdQ- The stable ground state \a) (|6)) is coupled 
to the excited state with a Rabi frequency 57i {^2) and 
a frequency uji (0^2) which is detuned from the excited 
state by Ai (A2). In the regime of weak laser power, 
the excited level is almost not populated and it can be 
adiabatically eliminated, leading to an effective photon- 
mediated interaction between the two ground state levels 
I a) and \b). By introducing the total angular momentum 

operators J+ = J2k=i l«)fc(^|fc' J- = J2k=i \b)k{a\k and 
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FIG. 4: Squeezing ^^ as a function of the cooperativity rj in 
the case of adiabatic (black circles) and optimal (red trian- 
gles) protocols with N = 30, 2k./ S ~ 10"^. The blue dashed 
line indicates the Heisenberg limit, reached in the large coop- 
erativity limit. 



"^z = iJ2k=i \a)k{a\k - |6)fe(&|fe)/2, and by further assum- 
ing the strength of the two Raman processes to be iden- 
tical, ftigl/Ai = ^29^/^2 = ^5*/A, after adiabatically 
eliminating also the cavity field, we obtain the following 
master equation for the density matrix |27| : 



p= -i[Heff,p] +Cp, 
with unitary part given by 

Heff 



X'JxJ 



(4) 



(5) 



where x = I^PlffP/'^^^ and S ^ ui — uio — uJab, and 
nonunitary part described by the Linbladian 



Cp = ^[2,P pJ- - J- .P p - pJ- J\ 



(6) 



where, from a microscopical derivation of the model, the 
most relevant contribution to the relaxation rate is 7 = 

Experimental implications — As discussed above, the 
Hamiltonian Eq. ([T]) is relevant e.g. for the experiment 
of Ref. [g . Here, squeezing of the collective spin of atoms 
in a cavity is used to improve the measurement precision 
of an atomic clock. With a realistic estimate of the pa- 
rameters [6| we have A = 780 nm '-^ 3 • 10^^ Hz, (5 ^ 27r • 3 
GHz, 7 - 27r • 5 MHz, 5 - 27r • 0.4 MHz, k - 27r • 1 MHz. 
The dominant part for the relaxation is thus proportional 
to the intensity of driving field x with a proportionality 
constant given by ^5/\g\'^ ^ 10^. Our estimate of relax- 
ation rate can be also expressed in terms of cooperativ- 
ity rj = 3^/(7^), leading to 7 = x(^)'^/(2t?7), where k 
is the decay rate of the cavity. In Fig. |4] the squeezing 
parameter is shown as a function of 77, for a system of 
size A'^ = 30. As clearly shown in the picture, the op- 
timized dynamics requires a cooperativity two orders of 
magnitude smaller than the one required by the (linear) 
adiabatic dynamics to reach the same level of squeezing. 



Conclusions and outlook — We have shown that opti- 
mal control can be used to speed up the dynamics for the 
production of squeezing, drastically reducing the effect of 
external dissipation. In particular we have demonstrated 
that optimized evolutions can occur on time scales or- 
ders of magnitude smaller than the corresponding adi- 
abatic evolutions, resulting in maximal squeezed states 
two orders of magnitude more robust with respect to the 
noise. The implementation of optimized protocols in spin 
squeezing experiments could therefore have a great im- 
pact in the field of quantum metrology. The very large co- 
operativity required here to achieve the Heisenberg limit 
is a consequence of a constraint intrinsic to the scheme 
of Ref. [J|, which aims at a pre-defined state (centered 
on top of the Bloch sphere as shown in Fig. [T]). Relax- 
ing this constraint on the position and orientation of the 
final state will allow to attain maximum squeezing al- 
ready with smaller cooperativity. The implementation 
of closed-loop optimal control strategies might result in 
additional improvement |25| . Finally application of the 
present methods demonstrated here to more complex 
spin squeezing schemes [5|, as well as adiabatic quantum 
computation in the presence of decohcrence, can also be 
envisioned. 
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Supplementary Material: Microscopical derivation of the model 

We consider a collection of N three level atoms with two stable ground states \a) and \b) and an excited state |e) in 
an optical cavity as done in. []|. The ground state \a) (|6)) is coupled to the excited state with a Rabi frequency fli 
{^^2) and a frequency uii (002), detuned from the excited state by Ai (A2). The laser frequencies are chosen such that 
their difference is equal to twice the ground state energy splitting, i.e. wi — W2 = '^'-^ab- With this choice all single 
atom transitions arc off-resonant and atoms can only be excited in pairs. The pairwise excitation of atoms is made 
possible by the presence of a quantized field inside the cavity coupling both states \a) and \b) to the excited state |e) 
with coupling constant ga and gt, see Fig.[T}D for the transition path involving a double Raman process. 
Assuming that all fields are propagating in the same direction, the Hamiltonian is: 
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H = UJqC^C + UJae^ee + ^ab^bb + 



ni 



-iiJit 



gac 



s.„, + ( ^ 



2 ^—iu>2t 



9bC 



h.c. 



(1) 



where S^j = X]a:=i l*}fc(j|fc are collective operators acting on the N atoms. 

In order to estimate the effect of the noise, we compute the equation of motion for the ground state coherence 
Sa6 = J+- The corresponding equation of motion takes the following form in a frame rotating at laser frequencies 
(H ^ H): 



n 



i'Eab = [^ah,H] = --—Y,),^ 



iv(i) 



2 ~eb ^5b*cte-'*'*I]W-ya.. -^6 ■ 3 



■5ace^'*sf 



m 



-^T^^a) - ho^ab + (Langevin force). 



where 70 is the dephasing rate of the ground state coherence, index 1(2) refers to left(right) Raman process in Fig.[Tl3, 
and the Langevin force (L.F.) term ensures the commutation relations of the operator. Analogously we can derive 
the equation of motion for the cavity field: 



-i[c,H] 



-^9be 



-^*I]W-^5:e-^''*Ei2i-^c+(L.F.), 



(2) 



with K being the decay rate of the cavity. 

In a regime of weak laser power, the excited level is almost not populated and it can be adiabatically eliminated. 

Therefore in such a regime where |fip/4 <C A^ + 7^, with 7 indicating the total decay rate of the excited state |e), we 

• 1(2) ' IC^) 

can set See ^ and Sea = ^he ^ ^- With these assumptions the equation of motion for the polarization becomes: 



iSafc = - 



"2 aI 



2 



-^Safc-.g,*cte-*'^*S6b 



9ace 
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— ^Sfcb - gac'e 



f~iSt 



Sat 



ni 1 + n/A2 






gbce'^'^ab 



^2^ 



gaCe Zjaa 



ijo^ab + (L-F.), 



X 



that is, separating the unitary part from the dissipative one 



iSab — [T,ab,H ] 
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(3) 



where H' is given by 
H' = 



4Ai 

ni 

2A^ 
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2 5,ce^^* I Sb 



(4) 



which corresponds to Eq. (Al) of Ref. (l|. The first line of previous equation represents ac-Stark shifts of the ground 
state: the classical part (proportional to f^'s) can be compensated making a change in the frequency of the fields; 
the part containing the quantum cavity field is much smaller than the first one and can be neglected in the limit 
g'^/SA << 1, see Ref. [ij. Notice that in this limit the part of the dissipation proportional to c'c can be also neglected. 
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FIG. 1: a) Left panel: Energy levels and couplings, b) Right panel: transition path involving a double Raman process. 



In a regime in which the lasers are sufficiently weak and we are not creating a significant number of photon excitations, 
the cavity field can be also adiabatically eliminated 
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dt 
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(5) 



where we are keeping only the Oth-order terms in the cavity field c. Notice that this equation corresponds to Eq. (A3) 
of Ref. [1| . Therefore inserting Eq. ^ in Eq. ([3]) and assuming the strength of the two Raman processes to be identical 
(fii<?*/Ai = f^2.9:/A2 = r!,g*/A), so that 
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1 — in/ 5 fig* 
S ^A 



(Eah + S(,a) 



we obtain 



il^ab = [Safe, H,ff] - z 70 + TTT + -TWT ) ^^b + i- ^^^, N (E^b + Sba) 



4A2 



4A2 



'A 2A2,5 



where we used E;,6 + 'Eaa '^ N and 



if, 



eff 






(6) 



(7) 



(8) 



Finally we can introduce the total angular momentum operators J_|- ~ Sat, J- = Sf,a and J^ = (Saa ^ '^bb)/'2, and 
the quantity x = jOplgp/fJA^ so that previous expressions become: 



iJ+ = [J+,Heff] 



M 70 + ^X ) J+ 



^^X^. + (L-F.) 



with the Langevin force term (L.F.) ensuring validity of the commutation relations for Jj-, and 

Heff = xJl- 



(9) 



(10) 



With a realistic estimate of the parameters [3 we have A = 780 nm ~ 3 • 10"'^'* Hz, 5 ~ 27r • 3 GHz, 7 '^ 27r • 5 
MHz, 5 -- 277 • 0.4 MHz, k - 277 • 1 MHz, so that we obtain for the imaginary part of Eq. dH) ^6/\g\'^ ^ 10^ and 
-fN/A ^ lO^'^N. The dominant part for the relaxation is thus proportional to the intensity of driving field x with 



a proportionality constant given by 7<5/|gp ~ 10^, so that wc can set in the master equation for the density matrix 
7 ~ xifyl^ l\g^ ^ as discussed in the main article. 
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